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Abstract 

We introduce new combinatorial (bijective) methods that enable us to 
compute the average value of three parameters of directed animals of a 
given area, including the site perimeter. Our results cover directed animals 
of any one-line source on the square lattice and its bounded variants, and 
we give counterparts for most of them in the triangular lattices. We thus 
prove conjectures by Conway and Le Borgne. The techniques used are 
based on Viennot's correspondence between directed animals and heaps 
of pieces (or elements of a partially commutative monoid). 

1 Introduction 

Let T be an oriented graph and S a nonempty finite set of vertices of T. A 
directed animal of source S on T is a finite set of vertices A that contains S and 
such that for every vertex v of A, there exists a vertex s of S and a path from s 
to v going only through vertices of A. The vertices of a directed animal A are 
called sites. The area of A, denoted by \ A\, is the number of sites of A. 

On Figure [T] are depicted directed animals on the three two-dimensional 
regular lattices: the square lattice, the triangular lattice, and the honeycomb 
lattice. The square lattice is the graph N 2 with the natural edges — > 
(i + l,j) and — > + 1) for all i,j > 0. The triangular lattice is simply 
obtained by adding the diagonal edges — > (i + l,j + 1). Traditionally, the 
lattices are depicted so that every edge points upwards. 






Figure 1: Single-source directed animals on the square, triangular and 
honeycomb lattices. All edges are directed upwards. 

Single-source directed animals constitute a subclass of animals (an animal 
on a non-oriented graph T is simply a finite connected set of vertices of F). 
While the enumeration of animals on any lattice is an open problem despite 
extensive research for decades, directed animals are fairly easier to enumerate. 
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As we will not deal with general animals in this paper, we will abusively use the 
term animal instead of directed animal. 

Single-source directed animals on the square and triangular lattices have 
been enumerated [H 12] • Specifically, let a(n) and a(n) be the number of 
animals of source {(0,0)} on the square and triangular lattice, respectively. 
The generating functions giving these numbers are: 

£«W = £^(7^-i); (i) 

n>l A \ " / 

Ya(n)t n = YtM = U- F ^-l). (2) 

In each case, the second sum goes over all single-source animals on the appro- 
priate lattice. 

Even then, much remains unclear. The enumeration of directed animals 
on the honeycomb lattice is an open problem; on the square and triangular 
lattices, comparatively very little is known when one tries to take into account 
parameters other than area. 



Today, two enumeration methods account for almost every known result on 
directed animals. One of them is the gas model technique, originally used by 
Dhar [B]. This technique was further developed and expanded by Bousquet- 
Melou [3] ; see also Q] for more recent work. 

The method used in this paper is the second one, based on a correspondence, 
due to Viennot [T2], between animals and other objects called heaps of dominoes. 
The basic idea is to replace each site of an animal by a 2x1 domino, so that 
each domino either lies on the ground or sits on one or two other dominoes 
(Figure^). 




Figure 2: A directed animal on the square lattice can be turned into a 
heap by replacing each site by a 2x1 domino. 



As we will see later, this method works for the triangular lattice as well. 
However, no simple model of heaps of dominoes has been found to correspond 
to animals on the honeycomb lattice. This may explain the lack of knowledge 
on the subject. 

The purpose of this paper is to study three other parameters of directed 
animals, introduced below, and illustrated in Figure U2 

• two sites of an animal on the square or triangular lattice are adjacent if 
they are of the form + and + We denote by j(A) the number 
of pairs of adjacent sites of A. 
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• a loop consists of two adjacent sites (i + l,j) and + 1), along with a 
third site at (i + l,j + 1). We denote by £(A) the number of loops of A. 

• a neighbour of an animal A of source S is a vertex v not in A, such that 
A U {v} is still a directed animal of source S. The number of neighbours 
of A is called the site perimeter of A and is denoted by p(A). 




( ) adjacent sites 

(^V) loop 

O neighbour 



Figure 3: A directed animal on the square lattice with two adjacent sites, 
a loop, and a neighbour marked. 



Taking, for instance, the site perimeter, we may consider the bivariate gener- 
ating function counting single-source animals according to both area and perime- 
ter on the square lattice: 

A*>(t,u)=Y, tlAluPiA) - 

A 

This generating function is not known, and is believed to be non-algebraic [9]. 
Instead, we will consider the generating function giving the total number of 
neighbours in the animals of a fixed area: 

A 

By dividing the total site perimeter of the animals of area n by the number of 
these animals, one gets the average site perimeter in animals of a fixed area. 
Alternatively, this generating function may be seen as counting single-source 
directed animals with a marked neighbour. 

This function, and the ones that similarly give the average number of adja- 
cent sites and loops, turned out to be easier to derive. Specifically, the value 
of the generating function counting the total number of loops of single-source 
animals on the square lattice was obtained by Bousquet-Melou using gas model 
methods [3]: 

As for the total site perimeter on the square lattice, it was the object of a 
conjecture by Conway in 1996 [5]: 

, .s mi 1 ( 2 1 - 3i + 2i 2 + t 3 - 3t i \ 

Le Borgne [10] also conjectured the value of similar generating functions count- 
ing the site perimeter of animals on bounded lattices. 
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Finally, we show in this paper that the total number of adjacent sites is given 

y^U-* fi- i-^+* 2 + 4 *;Y (5) 
Y zt(i+t)v v / r+t(i - 3*) 3 /v v; 

In Section HI we prove these three identities using combinatorial methods. Our 
formulae also instantiate in the case of animals of any one- line source, animals on 
the half- and bounded square lattices, and to a certain degree on the triangular 
lattices. All these lattices are defined in Section [3] 

Knowing, say, the total site perimeter of single-source animals of area n on 
the square lattice, we get their average perimeter by dividing by the number 
a(n) of these animals: 

a(n) 

y > \A\=n 

This quantity may thus be computed using ((T|) and Q . 

The numbers of adjacent sites and loops are handled similarly. From these 
generating functions, singularity analysis [7j yields estimates on these quantities 
as n tends to infinity: 

j(n) ~ -; £(n) ~ -; p(n) ~ — . 



The paper is organized as follows. In Section 2, we introduce in detail the 
notion of heaps of pieces and give several lemmas useful for animal enumeration. 
In Section 3, we enumerate directed animals of any source on several kinds of 
square and triangular lattices, according to area alone. In Section 4, we give a 
general method to derive the generating functions giving the average number of 
adjacent sites, number of loops, and site perimeter of directed animals on the 
square lattice; in Section 5, we give equivalents for most of these results in the 
triangular lattice. We illustrate our formulae with a few examples in Section 6. 



2 Heaps of pieces 
2.1 Definitions 

The notion of heaps of pieces, and its use to the enumeration of directed ani- 
mals, is due to Viennot, and this topic is covered in much more detail in |12j . 
Intuitively, a heap is a finite set of pieces. It is built by dropping successively 
the pieces at certain positions, chosen from a given set. When the position of 
two pieces overlap, the second piece falls on the first, like in Figure A more 
formal definition is given below. 

Definition 2.1.1. Let Q be a set and C a reflexive symmetric relation on Q. A 
heap of the model (Q,C) is a finite subset H of Q x N* verifying: 

1. if (<?, i) and (g', i) with q ^ q' are in H, then (g, q') is not in C; 

2. if (q, i) is in H and i > 1, then there exists (q 1 , i — 1) in H such that (q, q') 
is in C. 



4 



The relation C is called the concurrency relation, and two positions q and q' 
are concurrent if (q, q') is in C (in the above intuitive definition, this means that 
they overlap). 

The elements of a heap are called pieces. If (q,i) is a piece of a heap, q is 
called its position and i its height. The pieces of height 1 are called minimal; 
the set of the positions of the minimal pieces of a heap H is called the base of 
H, and denoted by b(H). A pyramid is a heap with only one minimal piece; 
a heap is trivial if all its pieces are minimal. As the height of the pieces of 
a trivial heap is meaningless, we will often confuse a trivial heap with the set 
of the positions of its pieces. In this way, a trivial heap is also a finite set of 
pairwise nonconcurrent positions. 

A piece (q, i) is said to sit on another piece of the form (q 1 , i — 1) if g and q' 
are concurrent. Thus, in Definition ^. 1.11 the second condition states that every 
non-minimal piece sits on at least one other piece. A finite set of pieces that 
verifies the first condition, but not necessarily the second is called a pre-heap. 

The reason we study heaps of pieces is to provide an alternate vision of 
directed animals on the square and triangular lattices as heaps of certain models, 
called heaps of dominoes. In these models, the families of heaps and pre-heaps 
defined below will correspond to directed animals on the square and triangular 
lattice, respectively. 

Definition 2.1.2. A heap (or pre-heap) H is strict if no piece sits on another 
piece with the same position. 

Definition 2.1.3. An inflated heap is a strict pre-heap H, such that for every 
piece (<7, i) in H such that i > 1, at least one of the following is in H: 

• either a piece (q',i — 1) such that q =/= q' and q and q' are concurrent; 

• or the piece (q, i — 2). 

With this definition, some non-minimal pieces of an inflated heap H may 
not sit on any other piece. We call these pieces unsupported; for reasons that 
will be clear when we link inflated heaps to animals on the triangular lattice, 
we denote by c(H) the number of unsupported pieces of H. 

Let H be a heap. A stack of pieces of H is a maximal set of pieces at 
the same position, sitting on each other (in other words, of the form (q, k) for 
* < k < i)- By applying the operations depicted in Figured] to each stack of a 
heap, we may link both strict and inflated heaps to general heaps. 

► 



Figure 4: A general heap may be obtained by replacing each piece of a 
strict heap by a stack of pieces; it can be turned into an inflated heap 
by inflating each stack. In an inflated heap, only the lowest piece of each 
stack is not unsupported. 
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2.2 Generating functions 

Now, we set some notations for the generating functions counting heaps of pieces. 
Let (Q,C) be a model of heaps and S a set of positions. 

Throughout this paper, we will use the same typographic pattern to denote 
generating functions counting heaps: the generating functions counting general 
heaps are denoted by calligraphic letters, while the ones counting strict heaps 
are denoted by standard capital letters. Moreover, the subscript S denotes 
heaps with base S, while the subscript [S] denotes heaps with base included in 
S. In this way, dehne the four generating functions: 

«s(t)= E <|ff| ; **(*)= E * m i 

b{H)=S b(H)=S 

H strict 

n [s] (t)= %](*)= E tim - 

b(H)CS b(H)CS 

H strict 

The generating functions counting heaps of base S and included in S are 
obviously linked by 

n [s] (t)= £ n T (t). (6) 

TCS 
T trivial 

To get the value ofHs{t) as a function of the Hm (t)'s, we may use the inclusion- 
exclusion principle: 

w s (t) = E(- 1 ) |s " T ' w m(*)- ( ? ) 

TCS 

Alternatively, define the neighbourhood of S, denoted by v(S) to be the set of 
positions concurrent to at least one position of S. Let H be a heap of base S. 
Then H consists of the pieces of S surmounted by a heap H', such that every 
minimal piece of H' sits on a piece of S. Thus, we have: 

H s (t)=t^H HS)] (t). (8) 

Identities similar to ^ and (J7J) also hold between generating functions counting 
strict heaps. 

We now translate the correspondences depicted in Figure 0] in terms of gen- 
erating functions. Replacing each piece by a nonempty stack is equivalent to 
the substitution 1 1— > j^-j- We thus get the following link: 

H s (t)=H s (-t-). (9) 



Or, equivalcntly: 



l-t 



//,(/) = //,( ^- ). (ID) 



Moreover, "inflating" a heap does not change its number of pieces, so that the 
generating functions counting general and inflated heaps are the same. Now, 
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let Hg(t,u) be the bivariate generating function counting inflated heaps, taking 
into account the number of unsupported pieces: 

b(H)=S 

As only the lowest piece of each stack of an inflated heap is not unsupported, 
this generating function is given by: 

m ,u)=Hs(j^). (11) 

We may now compute the generating function counting the total number of 
unsupported pieces of the inflated heaps with a given number of pieces: 

b(H)=S \ J \ / 

Using 0, this is equivalent to: 

c{H)t\ H \=t>H' s {t). (12) 

b(H)=S 

The case of heaps of base included in S is identical. 



2.3 The Inversion Lemma 

We say that a model (Q,C) is finite if the set of positions Q is finite. In this 
case, the heaps of any given base may be enumerated using the result below, 
due to Viennot [T2]. 

Let (Q,C) be a finite model and S a subset of Q. Define the following 
generating function, counting trivial heaps with positions included in S, giving 
a weight — t to each piece: 



I'i'l 



TCS 
T trivial 



As Q is finite, there are only finitely many trivial heaps included in S, and 
thus, this generating function is a polynomial, which is usually relatively easy 
to compute. 

Lemma 2.3.1 (Inversion Lemma). Let (Q,C) be a finite model, and let S be a 
subset of Q. The generating function 7i\s](t) exists and is given by 

J [Q] W 

In this way, we see that the generating function T~C\s] (t) is a quotient of two 
polynomials, and is therefore rational. 

In fact, the Inversion Lemma is valid as soon as the generating function 
T\q] (t) exists. It may exist even if Q is infinite, by giving the positions weights 
other than t, but we do not need it in this paper. 
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2.4 Monoid structure on heaps 

The main advantage of regarding directed animals as heaps of dominoes is the 
extra structure heaps provide. This structure is much more difficult to see on 
directed animals. To begin with, let H be a heap of a model (Q,C): then H is 
equipped with the following poset structure. 

Definition 2.4.1. Let if be a heap. Define the relation -< on the pieces of H 
by (q,i) -< (q',i') if and only if q and q' are concurrent and i < i' . Define the 
partial order < to be the reflexive transitive closure of -<. A piece x is said to 
be below a piece y (or y is above x) if we have x < y. 

This definition can be viewed more intuitively: let x be a piece of a heap H . 
Now, take the piece x and pull it upwards: the pieces that come up along are 
exactly the pieces above x. 

The minimal pieces of a heap H, as defined above, are actually those that 
are minimal with respect to <; those that are maximal for < are likewise called 
the maximal pieces of H. The set the maximal pieces of H is denoted by m(H). 
We say that two pieces of a heap are independent if neither is above the other. 

By turning a heap H upside down, letting all pieces fall in place, we obtain 
the dual heap H. This heap has the following property: a piece x of H is below 
another piece y if and only if x is above y in the heap H . In particular, the 
minimal pieces of H are the maximal pieces of H and vice versa. 

Let now H be a heap and q a position. The heap denoted by H ■ q is built by 
dropping a piece at position q on top of H. Formally, H-q is equal to HU{(q, i)} 
where i is the largest integer such that this is a heap. In this way, a heap may 
be built by dropping one piece at a time. This leads to the definition of the 
following operation on heaps. 

Definition 2.4.2. Let Hi and Hi be two heaps. The product Hi ■ H2 is built by 
dropping one by one all the pieces of H2 in any order compatible with <. 

The set of heaps of a model is a monoid under this law. Actually, this monoid 
is isomorphic to the partially commutative monoid on the alphabet Q with the 
concurrency relation C. An illustration of this product is given in Figure [5j 



I 



n — ► CD — ► □ 

□ □ □ cziz 



Figure 5: The product of two pyramids of dominoes, and the dual heap 
of the product. 

This structure is compatible with duality: let Hi and Hi be two heaps. 
Then the dual heap of the product H\ ■ Hi is Hi ■ H\ . 

2.5 Factorized heaps 

Definition 2.5.1. A factorized heap is a heap written as a product of two factors 
(Hi ■ Hi). If Hi ■ Hi = H, we say that (Hi ■ Hi) is a factorization of H. 
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Our goal is to give some results on factorized heaps which will be useful in 
the enumeration of animals. Specifically, if Hi is a fixed heap and S a trivial 
heap such that the base of H\ is a subset of S, we want to study factorized 
heaps of base S of the form (Hi ■ H 2 ). To do this, we need to figure out the 
base of a factorized heap (Hi ■ H2). 

Let H be a heap. Define the neighbourhood of H, denoted by v(H), to be 
the set of all positions concurrent to at least one piece of H (this definition 
extends the definition of Section [2~2l on trivial heaps). 

Lemma 2.5.2. Let Hi and H2 be two heaps. The base of the product Hi ■ H2 is 

b(Hi ■ H 2 ) = b(Hi) U (b(H 2 ) \ v(Hi)). 

Proof. Obviously, the minimal pieces of Hi are still minimal in Hi ■ H 2 ■ As for 
the minimal pieces of H 2l those that are minimal in Hi ■ H 2 are those not above 
a piece of Hi , thus the ones at positions not in the neighbourhood of Hi . □ 

Let now Hi be a heap, and let S be a trivial heap verifying b(Hi) C S. 
Let Tn^sit) be the generating function counting factorized heaps of the form 
(H x -H 2 ), such that b(Hi ■ H 2 ) = S. 

Lemma 2.5.3. Let Hi and H[ be two heaps sharing the same base and neigh- 
bourhood, and let S be a trivial heap such that b(Hi) C S. The following 
identity holds: 

^ H[ \s(t)=t^-^T Hl \ s (t). 

Proof. Consider the mapping (Hi ■ H 2 ) 1— > (H[ ■ H 2 ). As Hi and H[ share the 
same base and neighbourhood, Lemma [2 . 5 . 2l implics that it maps heaps of base S 
to heaps of base S . Moreover , we clearly have | H[ ■ H 2 \ = \Hi ■ H 2 \ + \H[\ — \Hi\. 
The lemma follows. □ 

We now consider the case of strict heaps. A direct analogue of Lemma [2. 5. 31 
does not hold, since the product of two strict heaps might not be strict. Instead, 
we use the following intermediate heaps. 

Definition 2.5.4. A factorized heap (Hi ■ H 2 ) is strict if Hi ■ H 2 is strict; it is 
almost-strict if Hi and H 2 are both strict. 

Let Hi be a strict heap, and let S be a trivial heap containing b(Hi). We 
denote by F Hl \g(t) the generating function counting strict factorized heaps of 
the form (Hi ■ H 2 ), with base included in S. We also denote by F^< s (t) its 
analogue counting almost-strict factorized heaps. 

Lemma 2.5.5. Let Hi be a strict heap, and S a trivial heap containing b(Hi). 
The generating functions counting strict and almost-strict factorized heaps of 
base S of the form (Hi ■ H 2 ) are linked as follows: 

F* Hlls (t) = (l + t)^ H ^F Hlls (t). 

Proof. Let (Hi ■ H 2 ) be an almost-strict factorized heap of base S. Then, in the 
underlying heap Hi ■ H 2 , the only pieces that might sit on another one at the 
same position are the minimal pieces of H 2 sitting on maximal pieces of Hi . 
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Let H' 2 be the heap formed by removing from H2 all minimal pieces at the 
same position as a maximal piece of H\. Then (Hi ■ H 2 ) is a strict factorized 
heap, and its base is still S. 

In this way, an almost-strict factorized heap is built from a strict factorized 
heap (Hi -H^) by possibly adding a piece under H' 2 at the same position as each 
maximal piece of H±. Each such piece "possibly added" accounts for a (1 + i) 
factor in the generating function F^, s , hence the formula. □ 

Lemma 2.5.6. Let Hi and H' x be two strict heaps sharing the same base and 
neighbourhood, and let S be a trivial heap such that b(Hi) C S. The following 
identities hold: 

^is(*) = * Wh|ffll ^i5(*); 

F H{\s{t) = q + t )|m(H»|-|m(ffi)| F Hi\s(t)- 

Proof. The proof of the first identity is similar to that of Lemma 12.5.31 as 
Hi and H[ have the same base and neighbourhood, Lemma f2 . 5 . 2 1 entails that 
the mapping (Hi ■ H2) 1— > (H[ ■ H2) is a bijection between the desired sets of 
almost-strict factorized heaps. The second identity is then a consequence of 
Lemma 12.5.51 □ 



2.6 Heaps marked with a set of pieces 

We are now concerned with heaps with a set of pieces marked. These objects 
intervene naturally in the enumeration of animals (for example, computing the 
average number of adjacent sites of animals of a given area is equivalent to 
enumerate animals with two adjacent sites marked). We define two mappings 
that link these heaps with factorized heaps. 

Definition 2.6.1. Let H be a heap marked with a set of pieces X. Define <d(H, X) 
to be the factorization (Hi ■ H2) of H such that Hi is the heap formed of all 
pieces below at least one piece of X (including the pieces of X). Symmetrically, 
define fl(H, X) to be the factorization (Hi ■ H2) of H such that H2 is the heap 
formed of all pieces above at least one piece of X. 

More intuitively, the factorization Q(H, X) is formed by taking the pieces of 
X and pulling them downwards, splitting the heap H into the desired two parts, 
while fl(H,X) is formed by pulling them upwards. Equivalently, VL(H,X) may 
be defined as the dual factorized heap of <d(H,X). 

We now consider the special case where the pieces of X are pairwise inde- 
pendent. 

Lemma 2.6.2. Both mappings 9 and f2 are bijections from the set of all heaps 
marked with a set of pairwise independent pieces to the set of all factorized 
heaps. Furthermore, let (Hi ■ H2) be a factorized heap. Its inverse image by 
is obtained by marking the set X of the maximal pieces of Hi ; its inverse image 
by f2 is obtained by marking the set of minimal pieces of i?2- 

Proof. Let (H,X) be a heap marked with a set of pairwise independent pieces, 
and let (H\-H^) be a factorization of H . We will prove that <d(H, X) = (Hi -H2) 
is equivalent to m(Hi) = X. 
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Assume that &(H, X) = (Hi ■ Hi). Let a; be a maximal piece of H\. As a; is 
in Hi, it is below one piece of X; as it is maximal, it is below no other piece of 
Hi than itself. Therefore, x must be in X. Let now y be a non-maximal piece 
of Hi: then it is below a maximal piece. As the elements of X are pairwise 
independent, y cannot be in X, so we have m(Hi) = X . 

Conversely, assume that m(Hi) = X. Every piece of Hi is below a maximal 
piece, and no piece of H2 can be below a piece of Hi in Hi ■ H^. Therefore, we 



In an enumerative context, these bijections are particularly useful, because 
the image of a marked heap if is a factorization of H. In particular, any heap 
has the same number of pieces as its image by 8 or (in fact, for every position 
p, the number of pieces at position p is also preserved). Thus, these bijections 
easily translate into identities between generating functions. 

Finally, we give a first practical application of the above results. Let S be a 
set of positions and q a position. Define the following two generating functions: 

• let Tvgj (t) be the generating function counting heaps of base included in 
S, marked with a piece at position q; 

• let Vjgi (t) be the generating function counting heaps of base included in S 
avoiding q, that is, such that no piece is concurrent to q (in other words, 
they are the heaps of base included in S of the model (Q \ v(q),C)). 

As usual, we use the same notations with standard capital letters to denote 
generating functions counting strict heaps. 

Proposition 2.6.3. The generating function counting heaps of base included in 
S marked with a piece at position q is given by: 



where H.{ p y(t) counts all pyramids of base {p}. 

Proof. Let H be a heap marked with a piece x at position q. We factorize the 
heap H by using the bijection O, pulling the piece x upwards: we thus have 
H = Hi ■ H-2, where H2 is the pyramid composed of all pieces of H above x. 

According to Lemma f2.5.2[ the base of H is b(Hi) if q is concurrent to a 
piece of Hi, and b(H{) U {q} otherwise. Thus, if q is in S then the base of Hi 
must simply be included in S; if q is not in S, then the base of Hi must be 
included in S and a piece of Hi must be concurrent to q. We therefore get the 
value of 

Let now Hy$*(t) be the generating function counting strict heaps of base 
included in S, with a piece at position q marked, such that the marked piece 
alone may sit on a piece at the same position. By factorizing such a heap in the 



have Q(H,X) = {Hi ■ H 2 ). 

The case of fi is easily obtained by duality. 



□ 




1 + 1 \ (h [s] (t) - Vjgj (t)) H {q} (t) otherwise, 
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same way, one gets an almost-strict factorized heap to which we may apply the 
same reasoning. Moreover, we have 



3 Directed animals and heaps of dominoes 

In this section, we will use the results on heaps of pieces to compute the gen- 
erating functions counting directed animals on various square and triangular 
lattices, according to area alone. To do this, we will define some models of 
heaps, called heaps of dominoes, and link them to directed animals. 

In Section [TJ we presented the full square and triangular lattices, which 
have the set of vertices N 2 (examples of animals on these lattices are shown in 
Figure [I}. We now define variants of these lattices, delimited by vertical lines. 

The half-lattice is delimited by one vertical line, while the rectangular lattices 
are constrained between two vertical lines. The cylindrical lattices are obtained 
by taking two vertical lines and gluing them together. 

The rectangular and cylindrical lattices are called bounded, as there is a 
bounded number of vertices on each horizontal line. The distance between the 
contraining vertical lines, which is always even in the cylindrical case, is called 
the width of the lattice (Figure EJ). 




Figure 6: Animals on the triangular half-lattice, and the square rectan- 
gular and cylindrical lattices of width 6. 

In addition, we restrain ourselves to animals with a one-line source, that is, 
a source which is a finite subset of a horizontal line. 

We denote by A(t) the generating function counting single-source animals on 
the full square lattice, and by D(t) the one counting single-source half-animals 
on the square lattice, that is, animals on the half-lattice with the source at the 
corner of the lattice. Moreover, let A2 m (t) be the generating function counting 
single-source animals on the cylindrical square lattice of width 2m, and D m (t) 
the one counting single-source half-animals on the rectangular square lattice of 
width m. 

In a way similar to heaps, we also use the notations Ait), T>(t), A2m(t) and 
T> m (t) to denote the generating functions counting the same animals on the 
various triangular lattices. 




(l + t)H§(t), 



as just one piece may be duplicated. 



This completes the proof. 



□ 
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3.1 Heaps of dominoes 

As stated before, we will compute these generating functions by regarding di- 
rected animals as heaps of dominoes. The model we consider depends on the 
kind of lattice: 

• the model Q = 1 corresponds to the full lattice; 

• the model Q = N corresponds to the half-lattice; 

• the model Q — [0, to — 1] corresponds to the rectangular lattice of width to; 

• the model Q = Z/2mZ corresponds to the cylindrical lattice of width 2m. 
In each case, two positions are concurrent if they are at distance at most 1: 



The strict heaps of these models are naturally in bijection with animals on the 
square lattices, while the inflated heaps are in bijection with animals on the 
triangular lattices. These bijections are illustrated in Figure [7] 

Let T be a square lattice, and A the associated triangular lattice. Let S be a 
one-line source: in terms of heaps of dominoes, S becomes a trivial heap, which 
we abusively also denote by S. Animals of source S on the lattices F and A 
become strict and inflated heaps of base S, respectively. Using the notations of 
Section |2~2"1 they are therefore counted respectively by Hs(t) and Tls(t). 

Definition 3.1.1. Let A be an animal of source S on the lattice A, and let 
be a site of A outside the source. We say that is only supported at the 

center if there are no sites of A at (i — or — 1); we denote by c(A) the 
number of sites of A only supported at the center. 

Let H be the heap of dominoes corresponding to A; a site v of A is only 
supported at the center if and only if the corresponding domino of H is unsup- 
ported. This enables us to use the results of Section l2~2l 

The heaps of dominoes obtained from animals also have an interesting prop- 
erty which we will exploit later. 

Definition 3.1.2. A (pre-)heap of dominoes H is aligned if the quantities q + i, 
for all pieces (q, i) of H, have the same parity. 

Let S be a one- line source. By construction, the trivial heap corresponding to 
S is always aligned. Moreover, the following result is easy to check by induction 
on the height of H: 

Lemma 3.1.3. Let H be a strict heap or inflated heap of dominoes such that its 
base b(H) is aligned. Then H is aligned. 

Thus, we conclude that the heaps of dominoes corresponding to directed 
animals on the square and triangular lattices are always aligned. 

In the next subsections, we give ways to compute the generating functions 
Hs(t) and Hs(t) counting animals of source S on the lattices F and A. We 
denote by a(n) and a(n) the number of directed animals of source S on these 
two lattices, so that: 



C = {(q,q')\\q-q , \<l}. 





n>0 



n>0 
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Figure 7: By replacing each site by a domino, animals on the triangular 
lattice become inflated heaps while those on the square lattice become 
strict heaps. In the animal at the top, two sites are only supported at the 
center. 



3.2 Bounded lattices 

Now, assume that T is a bounded lattice, cither rectangular or cylindrical. Then 
as the corresponding model (Q,C) of heaps of dominoes is finite, the Inversion 
Lemma (Lemma l2.3.1|) enables us to compute the generating functions counting 
the heaps of that model. 

Let T m {t) be the polynomial counting trivial heaps of the model [0, m — 1], 
each piece having weight —t. Likewise, let T2 m (t) be the polynomial counting 
trivial heaps of the model Z/2mZ. The following proposition gives an inductive 
formula to compute these polynomials. 

Proposition 3.2.1. The polynomials counting trivial heaps of dominoes on the 
finite models with each piece having weight —t are given by: 

if m < 0; 
otherwise; 

Actually, these polynomials are well-known: they are closely related to the 
Fibonacci (or Chebyshev) polynomials. 

Sketch of the proof. Let T be a trivial heap of the model [0, m — 1]: then T is 
counted by T m _i(t) if m — 1 is not in T, and by — iT m _2(i) otherwise. The 
cylindrical case is handled similarly. □ 

Using this result and the Inversion Lemma, we are able to derive the values 
of the generating functions counting general (or equivalently, inflated) heaps of 
dominoes of any base, hence animals on the triangular bounded lattices. For 
example, we give below the generating functions counting single-source half- 
animals and animals. 



T m {t) = 1 

T m (t) = T m _i(t) — tT m ^2(t) 

T2m(t) = T2m-l(t) ~ tT2 m —3(t)- 
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Proposition 3.2.2. The generating functions for single-source animals on the 
cylindrical triangular lattice and half-animals on the rectangular triangular lat- 
tice are: 

A 2m (t) = %^ - 1; (13) 

The generating functions counting the corresponding animals on the square 
lattices are found by replacing t by j-fj- Using this result, we are able to derive 
the asymptotic behaviour of the number of animals of area n as n tends to 
infinity. 

Lemma 3.2.3. For all m > 0, The polynomials T m {t) and T 2m (t) have only real, 
simple roots. Let p m and ai m be their smallest respective roots. The sequences 
{pm) and (<72m) are decreasing and tend to 1/4 as m tends to infinity; moreover, 
we have aim < Pim f° r all m. 

Proof. First, we check by induction using the definitions of T m (t) and T 2m (t) 
the following identities: 

1 \ sin[(m + 2)0] 



1 2m 



4 cos 2 9 J (2cos0) m + 1 sin(9' 
1 \ 2cos(2m0) 



4cos 2 0/ (2cos0) 



2m ' 



We also check that the polynomials T m (t) and T 2m (£) have degree [ttJ and to, 
respectively. This proves that all roots of both polynomials are real and simple; 
moreover, we have: 



4 rOS 2 7T ' 

^ t -° fe m+2 



C2m 



■1 COS 2 ^ ' 



The lemma follows from these values. □ 

Proposition 3.2.4. Assume that T is a bounded square lattice and let A be the 
associated triangular lattice. Let S be a one-line source. As n tends to infinity, 
the number of animals of source S and area n on the lattices T and A are 
asymptotically of the form: 

a(n) ~ A// 1 ; a(n) — \p, n . 

Moreover, while the constants A and A depend on the source, p depends only 
on the lattice, and p is equal to p + l. 

In particular, if S and S' are two one-line sources, the number of animals of 
area n with source S and S' on any bounded lattice differ only by a multiplicative 
constant. 



15 



Proof. Assume that T is the rectangular lattice of width m. We derive the gen- 
erating function 7iis]{t) using the Inversion Lemma: its numerator is a product 
of factors Tfc 4 (i), for some ki < to, and its denominator is T m {t). We then use 
singularity analysis: Lemma 13 . 2 .31 entails that the smallest pole of Tt[s](t) is p m , 
and that this pole is simple. We thus get the estimate for a„, with p — l/p m - 

The case of the cylindrical model of width 2m is identical, with the denom- 
inator being T2 m (t). As all pk for k < 2m are greater than <72 m according to 
Lemma 13.2.31 we find /2 = l/c2 m . 

The generating function H[s](t) is obtained by substituting t by j^t- Thus, 
the smallest pole 1/p of H[s](t) verifies 

1/M = 1 
1 + p, 

which yields p, = fj, + 1. As l/ft is a simple pole of Hrsi(t), this pole is also 
simple. 

The case of animals with source S is handled with the inclusion-exclusion 
principle. □ 



3.3 Unbounded lattices 

Now, we deal with the case where T is an unbounded lattice; as usual, we 
denote by A the associated triangular lattice. The generating functions counting 
animals and half-animals on the square and triangular unbounded lattices are 
given below. Various ways to derive them can be found in O [8j [2] ■ 

Proposition 3.3.1. The generating functions counting half-animals and animals 
on the unbounded square and triangular lattices are given by: 




- v/(l + t)(l-3f) 



(15) 



Vi- 


-4t 


2t 




ft 




-3t 


1 




Vi- 


At 



2t 

-1; (16) 



-i ; (17) 



-1 . (18) 



Of course, as usual, the generating function counting animals on the square 
and triangular lattices are linked by: 



V{t) = D 
A{t) = A 



t 



1 - 1 
t 

1 - 1 



Moreover, we have, as to tends to infinity, the following limits in the space of 
formal power series: 

V m (t) ->2>(t); 
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If S is a one-line source, we will now compute the generating functions 
counting animals of source S on the lattices T and A. To do so, we see these 
animals as strict and inflated heaps of base S on the corresponding model Q of 
heaps of dominoes. 

In the next propositions, we give only the generating function counting gen- 
eral, or inflated, heaps of base S. To get the one counting strict heaps, it suffices 
to replace t by YTt- Thus, the generating functions V(t) and A{t) are replaced 
by their non-calligraphic counterparts. 

In the model Q = Z, for any position q, the pyramids of base {q} are counted 
by the generating function A{t). In the model Q = N, this is not the case, as 
shown by the following result. 

Lemma 3.3.2. Let Q = N be the model corresponding to the half-lattice, and let 
q be a position. The generating function counting pyramids of base q is given 
by 

^ (t) " i-v{t) v{ty 

Regarded as directed animals, these pyramids are more commonly called 
"directed animals with left width at most q" . A proof of the above formula is 
found in [2]. 

We call a trivial heap compact if it is of the form {q, q + 2, . . . , q + 2k — 2}, 
with k > 1. The heaps with a compact base are enumerated by the following 
proposition. 

Lemma 3.3.3. Let Q — Z and let Ck be a compact trivial heap with k pieces. 
The generating function counting heaps with base Ck is given by: 

Hc h {t) = A{t)V{t) k - 1 . 

Let now Q = N and let C q ^ be the compact trivial heap {q, q + 2, ■ ■ ■ , q+2k — 2}. 
The generating function counting heaps with base C q .k is given by: 

By summing the first identity over all k > 1, we get this surprisingly simple 
result, also found in [S]. 

Proposition 3.3.4. The number of directed animals with any compact source on 

the square lattice with area n is 3™ _1 . On the triangular lattice, this number is 
4«-l_ 

Proof of Lemma 13.3.31 To prove this result, we use the construction depicted in 
Figure [8J due to Betrema and Penaud [2] . 

In this way, a heap with base C q ^ can be seen as the product of k — 1 half- 
pyramids and one pyramid with base {q}. This yields the announced generating 
functions. □ 

Finally, we are able to compute the generating function counting heaps of 
base S, for any aligned trivial heap S. 
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Proposition 3.3.5. Let Q be either infinite model of heaps of dominoes (Q = Z 
or Q = N), and let S be an aligned trivial heap. Let C be the smallest compact 
trivial heap such that S C C. The generating function counting heaps with 
base 5* is given by: 

Hs(t) = Tlv % v J°/ t] H c (t) 

where T\ v tc)^v(S)](t) counts all trivial heaps included in v(C) \ v(S), each piece 
having weight —t. 

In the case of the full lattice, this result is found in an equivalent form in pQ. 

Proof. First, define the finite model Q m in the following manner: if Q = Z, let 
Q m be the model Z/2mZ; if Q = N, let Q m = [0, to — 1]. In either case, m is 
assumed to be large enough so that S and C are still trivial heaps of the model 

{Qmi C ) • 

Using identity (|8|) and the Inversion Lemma, we compute the generating 
functions counting heaps of this model with base <S* and C: 

T~Lc,m(t) 

As S is included in C, v(S) is also included in v(C), and therefore, we can write: 

Q m x v[S) = (Q m \ v(C)) U (v(C) \ v(S)). 

As C is a compact source, v(C) is a segment of integers, so that the only positions 
of v{C) concurrent to a positions of Q m \ v(C) are its endpoints. As C is the 
smallest compact source containing S, both endpoints of v(C) are also in v(S). 
Therefore, no position of Q m \ v(C) is concurrent to a position of v{C) \ v(S). 

Thus, we see that choosing a trivial heap outside v(S) is equivalent to 
choosing independently two trivial heaps, one outside v(C) and the other in 
v{C) n v(S). In terms of generating functions, we have: 

^mNtfSilW = T[Q m ^ v ( C )](t)7[ v (C)^v(S)](t)- 

Therefore, we get: 

We conclude by letting m tend to infinity. □ 



= * 1[,(S)],m(')=f 



^t 1 'rt^c ,mW-t' 7f TTn • 

7 [<2m]W 
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Similarly to the case of bounded lattices, we get the asymptotic behaviour 
of the numbers a(n) and a(n) of directed animals of source S and area n on the 
lattices T and A. 

Proposition 3.3.6. Assume that T is the full lattice, and let S be a one-line 
source. As n tends to infinity, we have: 

3™ - 4™ 

a(n) ~ A— a(n) ~ A— p= 
V n v n 

where A and A are constants depending on the source S. 

Assume now that T is the half-lattice. The number of animals of area n is 
equivalent to: 

3™ - 4" 

with A and A again depending on the source S. 

Again, if S and S' are two one-line sources of a lattice, the number of animals 
of source S and S' of area n asymptotically differ only by a multiplicative 
constant. 

Proof. Once again, we use singularity analysis. The smallest pole of the gener- 
ating functions A(t) and T>(t) is t = 1/4, and we have as t tends to 1/4: 

A{t) = (1 - At)- 1 / 2 + o((l - At)- 1 ' 2 )- 
V(t) = 1 - 2y/T^M + o(V i - it) . 

Now, assume that Q = Z and let Cfc be a compact source. According to 
Lemma l3.3.3[ we have: 

W Cfc W = (l-4i)- 1/2 +o((l-4r 1/2 ). 
If Q = N, then let C 9i fe be a compact source. We have: 



K c ,, t (t) = q + 1 - (g + l)(2fc + ^VI - 4< + o(VT^ 4t). 

Hence, we get the announced estimates if 5 is a compact source. 

Let S be any one-line source, and let C be the smallest compact source 
containing S. We use Proposition I3.3.E>1 to link the generating functions Tis(t) 
andTtcW- Thus, t\ c \-\ s \n s (t) /H c {t) is a polynomial counting trivial heaps 
of a finite model, which according to Lemma 13.2.31 has all roots larger than 
1/4. Hence, this polynomial accounts only for a multiplicative constant in the 
behaviour of the coefficients of Tts(t). 

As in the proof of Proposition ^. 2. 41 the case of the square lattice is handled 
by performing the substitution 1 1— > j^~t ■ ^— ' 

4 Average number of adjacent pieces, number of loops, 
and site perimeter on the square lattice 

4.1 Notations and results 

In this section, T is a square lattice (full, half, cylindrical or rectangular), and 

5 is a one-line source of T. In Section [31 we have computed the generating 
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functions Hs(t) and H^]{t) counting directed animals on T of source S and 
with a source included in S. 

In Section 1, we have defined the number j(A) of pairs of adjacent sites and 
£(A) of loops of an animal A. Now, define the generating function, counting the 
total number of adjacent sites of animals of a given area: 

Js(t) = '£j(A)tW 

A 

where the sum goes over all animals of source 5 on T. In the same way, define 
the generating function Ls(t) counting the total number of loops; also define 
the analogous generating functions J[s](t) and L[s](t) dealing with animals of 
source included in S. 

Definition 4.1.1. Let A be an animal on T with a source included in S. A 
neighbour of A is a vertex v of T not in A such that A U {v} is still a directed 
animal with a source included in S. 

Moreover, assume that the graph T is embedded in a larger graph V . An 
internal neighbour of A is a neighbour of A seen as an animal on T. An external 
neighbour of A is a neighbour of A seen as an animal on V. 

For the purpose of this definition, we regard the half-lattice and the rectan- 
gular lattices as embedded in the full lattice. The full and cylindrical lattices 
are not naturally embedded in any other graph than themselves, so in these 
lattices internal and external neighbours are identical. 

We denote by Pi(A) the number of internal neighbours of A (or internal site 
perimeter) and by p e {A) its number of external neighbours (or external site 
perimeter). Thus, we define the generating functions Ps(t), -PJ(i), P[ S ](t) an d 
Pj^j (t) in the same way as above. 

To compute these generating functions, we again use the correspondence 
between animals and heaps of dominoes. We denote by (Q,C) the model of 
heaps of dominoes associated to the lattice T; we also denote by S the aligned 
trivial heap associated to the source S. 

In this section, we express all the above defined generating functions in 
terms of other generating functions that can all be computed using the results 
of Section [3] Namely: 

• let Es (t) be the generating function counting strict heaps of base S marked 
with a piece at the edge of the model, i.e. at a position q such that either 
q — 1 or q + 1 is not in Q; 

• let Us(t) be the generating function counting strict heaps of base 5* marked 
with a piece at a position q such that q + 2 is not in Q; 

• if q — 2 is in S but q is not, let Wg(t) be the generating function counting 
strict heaps of base S U {<?}; 

• let W s (t) be the sum of W q s {t) for all q in (S + 2) \ S. 

Remark. As the full and cylindrical models have no edge, the generating func- 
tion E$(t) is always zero in these models. Moreover, in all models but the 
rectangular, the position q + 2 is always in Q if q is, so the generating function 
Us(t) is also zero. 
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In the same way, define the generating functions (t) and U[g] (t) ■ Finally, 
define: 

• if q — 2 is in S but q is not, is the generating function counting 
strict heaps of base included in S U {q} containing q — 2 and q\ 

• W [s] (t) is the sum of WL(t) for all q in (S + 2) \ S. 

The results of Section[3]enable us to compute the generating functions Ws(t) 
and W[g](t). Moreover, as there are at most three positions q such that q — 1, 
q + 1 or q + 2 is not in Q, we can compute the generating functions Eg(t), Ug(t) 
and their counterparts using Proposition 1 2 . 6 . 3l 

The following theorem, proved in Sections 14.21 and |4.3[ thus gives the value 
of the generating functions counting the total number of adjacent sites, loops 
and site perimeter of directed animals on T. 

Theorem 4.1.2. The generating functions counting the total number of adjacent 
pieces, loops and site perimeters of the animals with source included in S are 
given by: 



t*H' [S] (t)-tU lS] (t)-W [S] (t) 

J is]{t) = j— - t ; (iyj 

L [s] (i)=t(i + t)J [s] (i); (20) 

P( s] (t) - \S\H [S] (t) + tH' [s] {t) - J [s] (t); (21) 
^S](*)=-P[S](*) -%](*)• (22) 
Moreover, the analogous generating functions for animals of source 5* are: 
Js {t) = t 2 H' s (t) - tU s (t)+j(S)H s (t) - W s {t) , 

1 "I - ^ 

L s (t)=t(l+t)Js(t); (24) 
P|(t) = |S|ff 5 (t)+t^(t)-^(*); (25) 

Pg(t) = Pl(t) - Eg(t), (26) 



where ^(S 1 ) denotes the number of pairs of adjacent pieces in the trivial heap S. 

In Section [SJ we give applications of this theorem to several specific lattices 
and sources. 

4.2 Links between adjacent sites, loops and site perimeter 

Our first step to prove Theorem 14.1.21 is to show that the total number of 
adjacent sites, loops and both site perimeters of directed animals are linked. 

We begin by showing a link between adjacent sites and the site perimeters. 
If A is an animal, let e(A) be the number of sites of A at the edge of the lattice, 
i.e. the number of sites with outgoing degree 1 in T. Thus, we have: 

E s (t) = '£e(A)tW, 

A 

where the sum goes over all animals A of source S. 
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Lemma 4.2.1. Let A be an animal with a source included in S on L. The 
external and internal site perimeters of A verify: 

Pe (A) = \S\ + \A\-j(A): 
Pi (A)= Pe (A)-e(A). 

By summing these identities over all animals with a source included in S 
and of source S, we get the identities (|2"T|) . ((HJ) and (ggj) . 

Proof. When dealing with the external site perimeter, the lattice L is embedded 
in a lattice V which is either the full lattice or a cylindrical lattice. Let Z be 
the number of pairs of vertices (v, w) such that v is a site of A and w is a child 
of v (i.e., v — > w is an edge of T), whether in A or not. As every vertex has 
outgoing degree 2, we have 

Z = 2\A\. 

Now, as A is a directed animal, a child of a site of A is either a site of A or 
a neighbour of A. The only sites and neighbours not counted are the ones in 
5 1 ; moreover, two sites have a child in common if and only if they are adjacent. 
Hence: 

Z = \A\+p e (A)+j(A)-\S\ 

which yields the announced formula for p e (A). 

If F is either the half-lattice or a rectangular lattice, then each site on the 
edge of the lattice has one external neighbour not in L. Thus, we have 

p e (A)= Pi (A) + e(A). □ 

The next result links the number of adjacent sites to the number of loops, 
and yields the identities ([2U]) and ■ 

Lemma 4.2.2. The generating functions counting the total number of adjacent 
sites and loops verify: 

L s {t)=t(l + t)J s (t); 
L [s] (t)=t(l + t)J [s] (t). 

Proof. Let A be an animal with two adjacent sites marked. We regard A as a 
strict heap of the model (Q,C) marked with two adjacent pieces. By pulling 
these two pieces downwards, we get a factorized heap (Hi ■ H2), with base 5* 
and such that H\ has two maximal, adjacent pieces. 

We do the same to an animal marked with a loop, and we get a factorized 
heap (H[ ■ H2) with base S such that H[ has one maximal piece sitting on two 
adjacent pieces. 

Placing a piece on two adjacent pieces does not alter the base or neighbour- 
hood of Hi , so the first identity is a direct application of Lemma 12.5.61 The 
second easily follows by summing over all sources T C S. □ 
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4.3 Average number of adjacent sites 



With the results of Section 14.21 all the identities of Theorem 14.1.21 are con- 
sequences of the equations (JH]) and (|23p . dealing with the average number of 
adjacent pieces. First, we concentrate on (|19j) . considering animals with a source 
included in S. We will later derive (|23[) using an inclusion-exclusion argument. 

To count animals marked with two adjacent sites, we regard them as heaps 
marked with two adjacent pieces. We then pull these two pieces downwards, cre- 
ating a factorized strict heap (Hi ■ H2), such that Hi has two adjacent maximal 
pieces. 

In order to handle these factorized heaps, we consider almost-strict factorized 
heaps rather than strict. Define the following sets of factorized heaps: 

• J* the set of almost-strict factorized heaps (Hi ■ H2) such that Hi has 
exactly two maximal pieces, which are adjacent; 

• M* the set of almost-strict factorized heaps (Hi ■ H2) such that H\ has a 
unique maximal piece at a position q such that q + 2 is in Q; 

• I 2 * (resp. I 3 *) the set of factorized heaps (H\-H2) such that Hi has exactly 
two maximal pieces, at some positions q and q + 2 (resp. q and q + 3). 

Moreover, let J* S ](t), M^(t), lfgt(t) and I?g](t) be the generating functions 
counting factorized heaps of these sets with base included in S. 

Also, if q — 2 is in S but q is not, let W^Jt) be the generating function 
counting almost-strict factorized heaps (Hi ■ H2), with base included in SU {q} 7 
such that Hi is reduced to a single piece at a position q — 2 and H2 has at least 
one minimal piece at position q. Let W*<^ (t) be the sum of W?Z (t) for all q in 
(S + 2)\S 

Lemma [2.5.51 gives the link between almost-strict and strict factorized heaps: 

J{ s] (t) = (l + t) 2 J [S] (t); (27) 
M[S] (*) = (1 + 1) (tH' [s] (t) - U [s] (t)) ; (28) 
W [ * s] (t) = (l + t)W [s] (t). (29) 

Lemma 4.3.1. The following identity holds: 

Ifs](t) = -fa (*)+*/&(*)• 
Proof. The core of this lemma is the bijection 

$ : I 2 * \ J* -> I 3 * 

illustrated in Figured! Let (Hi ■ H2) be a factorized heap in I 2 * \ J*. Therefore, 
the heap Hi has two maximal pieces at some positions q and q + 2, and these 
pieces are not at the same height. Let x be the lower of the two, and y the 
higher. 

The heap $(Hi ■ H2) is obtained by removing the piece y. This piece cannot 
be minimal, and as Hi is strict and a; is a maximal piece, y must sit on a piece 
z on the opposing side of x. Thus, • H2) is a heap of I 3 *. 

Conversely, let (Hi ■ H2) be a factorized heap in I 3 *. The heap Hi has 
two maximal pieces at some positions q and q + 3. As Hi is strict, it is aligned 
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I 2 * \ J* 

Figure 9: The bijection after removing the piece y, a maximal piece z 
is uncovered at a position at a distance 3 from x. 

according to Lemma f3. 1.31 Therefore, the maximal pieces' heights have opposite 
parity and thus cannot be equal. Let x be the lower piece and z the higher: 
then $ _1 (Hi ■ H 2 ) is obtained by adding a piece y on z at position q+1 or q + 2. 
In addition, let (Hi ■ H 2 ) be a factorized heap of I 2 * \ J* and let (H[ ■ H 2 ) = 
• #2)- As the heaps Hi and have same base and neighbourhood, 
Lemma \2 . 5 . 21 entails that the bijection <!> preserves the base. Moreover, it satis- 
fies: 

\<S>(Hi ■ H 2 )\ = \H X ■ H 2 \ - I. 
This yields the announced identity on generating functions. □ 
Lemma 4.3.2. The following identity holds: 

I[s) (*) = tlfs) (*) + tM*s] (*) " ^[S] (*)■ 
Using this lemma and the previous one, we derive (fT§|) using |2T|) . (|28p and 

Proof. Again, we prove this lemma using a bijection 

* : I 2 * -> I 3 * U M*. 

This bijection is illustrated in Figure fTOl Let (.Hi • ^2) be a factorized heap of 
I 2 *. The heap Hi has two maximal pieces x and y, at positions g and q + 2 
respectively. 

Again, the factorized heap ^(iJi ■ H 2 ) is obtained by removing the piece y. 
Thus, the heap Hi \ {y} may either have two or one maximal pieces. As Hi is 
strict, in the former case, the new maximal piece z must be at position q + 3, 
so we obtain a heap of I 3 * or M* respectively 




I 3 * U M 



Figure 10: The bijection ^: after removing the piece y, either a maximal 
piece z is uncovered at position q + 3 or not. 

The inverse mapping ^E' _1 (_ffi • H 2 ) is obtained by adding back a piece y on 
top of Hi at position q + 2, where q is the position of x. The bijection ^ again 
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satisfies: 

\^{H X -H 2 )\ = \H 1 -H 2 \-\. 

However, we cannot directly restrict the bijection 'F to heaps of base included 
in S because it does not preserve the base. There are two cases when the base 
of ^>(Hi ■ H 2 ) is different from that of H\ ■ H 2 , illustrated in Figure ITTl 

1. if the piece y is a minimal piece of Hi and H 2 has no minimal piece at 
position q + 2, then q + 2 is in b{H x ■ H 2 ) but not in • H 2 )) ; 

2. if there is a minimal piece z of H 2 with position q + 'i and there is no piece 
of Hi concurrent to it, then q+3 is in b(Hi-H 2 ) but not in b( i l>~ 1 (Hi-H 2 )}. 

Let (t) be the generating function counting the factorized heaps (Hi-H 2 ) 
with base included in S and verifying condition 1 above, such that q + 2 is not in 
S, and let q' — q + 2; let Y? s At) be the one counting those verifying condition 2, 
such that q + 3 is not in S, and let q' = q + 3. The bijection ^ thus entails that 

41 (*) + x k (*) = *( J is] (*) + M k (*) + 3% (*)) ■ 



112 ;-x y ^ 2 m / 




Figure 11: At the top, the two cases where the bijection does not 
preserve bases, counted by X? 3 -,(t) and Yjm(t). A heap counted by Xfm(t) 
is either counted by W^(t) or by tYrg](f), depending on whether Hi is 
reduced to the pieces x and y or not. 

To simplify this equation, we perform the operations depicted in Figure ITTl 
Let (Hi ■ H 2 ) be a factorized heap counted by X\s\(t); let Hi = H[ ■ x ■ y, and 
H 2 = yH 2 : 

• if H[ is empty, then the factorized heap (x ■ H' 2 ) is counted by W,* s Jt); 
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• if not, then as y is a minimal piece of H±, no piece of iJ{ is concurrent to 
q'\ therefore, as H± is strict, H[ has a unique maximal piece x' , at position 
q' — 3 (Figure fTTl bottom right). Therefore, the factorized heap (H[ ■ H 1 ^) 
is counted by Y* s ^(t). 

Finally, we find: 

X(s ] (t) = W^ s] (t) + tY^ ] (t). 
This completes the proof. □ 

Having proved the identity (|19|) . we now derive (|23|) using the inclusion- 
exclusion principle. This principle states that if Z^(t) and Z$(t) are generating 
functions verifying 

Z [s] (t)=J2 Z T(t), 

TCS 

then Zg(t) is given by: 

TCS 

The generating functions Htg](t), J[s](t) an d ^[S] W n t this description. To 
deal with W[g](t), we need the lemma below. 

Lemma 4.3.3. The following identity holds. 

£ (-1) |S " T| ^[T]W = W s (t) -j(S)H s (t). 

TCS 

With this lemma and the inclusion-exclusion principle applied to J[s](t), 
H' [s] (t) and ?7[s](t), the identity (J23J) follows from [jig]). We have thus finished 
proving Theorem 14. 1.21 

Proof. Consider the generating function 

E ^(t). 

TCS 

Recall that Wx(t) is the sum of all W T (t), for all q such that q — 2 is in T but 
q is not. We further distinguish whether q is in 5*: 

E wtw = £ ( £ + E f E *l w) ■ 

TCS TCS\ 9 £S / TCS\ gCS / 

q-2CT 9-2GT 
9 £T 

In the first term of the r.h.s., we recognize the generating function W[s](t). We 
rewrite the second term using the fact that W T (t) = H Tu ^ q y(t) and by posing 
T' = TU{q}: 

E w T (t) = w [s] (t) + e ( E h t- (*) J ; 

TCS T'CSX gG T' / 

g-2CT' 

= Wfa(t) + E 3{T')H TI {t). 

T'CS 
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We rewrite this as: 

W [s] (t)=Y J {W T (t)- ] (T)H T {t)). 

TCS 

Thus, we may apply the inclusion-exclusion principle on T4 / [sj(t), which yields 
the announced formula. □ 



5 Triangular lattices 

Here, we state some results on animals on the triangular lattices, similar to those 
of Section 0] to a certain extent. We consider the triangular lattice A associated 
to r, and we are interested in the average number of adjacent sites, loops and 
neighbours of animals of a given area. 

Adjacent sites and loops are defined in Section [T] Note that the definition 
of loops we use is somewhat unusual; it is, for instance, different from the one 
found in [3] . Like on the square lattice, we now introduce several notions of site 
perimeter. 

5.1 Site perimeter 

Let A be an animal with a source included in S. We denote by Pi(A) and p e (A) 
the internal and external site perimeters of A, defined in a manner identical to 
Definition 14.1.11 

Unlike on the square lattice, the generating functions giving the average 
value of the site perimeters are believed to be non-algebraic, and we have not 
been able to compute them. Instead of studying them directly, we consider the 
two quantities below. Let c(A) be the number of sites of A only supported at 
the center (see Definition ^. Lip . 

Definition 5.1.1. Let A be an animal with a source included in S. Define the 
number p_ (A) to be the number of external neighbours v of A, such that v is 
not only supported at the center in the animal A\J {v}. Moreover, define: 

p + (A) = 2\S\ + \A\+c(A)-l(A). (30) 

Let e(A) be the number of sites of A at the edge of the lattice A (that is, 
vertices having outgoing degree two). 

Lemma 5.1.2. Let A be an animal with a source included in S. The external 
and internal site perimeters of A verify: 

p.(A) = \S\ + \A\+c(A)-j(A); 
P-(A) < p e (A) < p+(A); 
Pl (A)=Pe(A)~e(A). 

Proof. The proof is very similar to that of Lemma 14.2.11 on the square lattice. 
First, we consider the external site perimeter and therefore assume that the 
animals live in an extended lattice A' where every vertex has outgoing degree 
three. 
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Let us begin with p_(A). Let Z be the number of pairs (v, w) of vertices of 
A', such that v — is a site of A and w is of the form + 1) or (i + 
Obviously, we have 

Z = 2\A\. 

Likewise, if (v, w) is such a pair, then w is either a site of A or a neighbour of 
A. The only sites and neighbours not counted are the vertices of S and the sites 
and neighbours only supported at the center. Moreover, a vertex is counted 
twice every time two sites of A are adjacent. Therefore, we have: 

Z = \A\+p_(A)-\S\-c(A)+j(A). 

The inequality p-(a) < p e {A) is obvious by definition. 

Now, let us prove the upper bound on p e (A). Let Z' be the number of pairs 
(v, w) of vertices of A', such that v is a site of A and w is a child of v. As every 
vertex has outgoing degree three, we have: 

Z' = 3\A\. 

For each such pair, w is either a site or a neighbour of A, and the only sites and 
neighbours not counted are the vertices of S. Moreover, as shown in Figure \12\ 
a vertex w is counted twice at least every time a site v of A is supported on the 
left or on the right. 



w w 




Figure 12: Each time a site v is supported on the left or on the right, a 
vertex w is the child of two sites. 

Let v be a site of A. There are three possible cases: 

• the site v is neither supported on the left nor on the right, i.e. is either in 
S or only supported at the center; 

• the site v is supported both on the left and right, i.e. is a loop; 

• the site v is either supported on the left or on the right, but not both. 
Thus, we get: 

Z' > \A\ + Pe (A) - \S\ + 21(A) + (\A\ - \S\ - c(A) - 1(A)) 
which, by the definition ()30|) . boils down to 

p e (A) <p+(A). 

Finally, the link between external and internal site perimeter is due to the 
fact that every site of A with outgoing degree two in A has one less neighbour 
in A than in A'. □ 
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5.2 Average number of adjacent sites, number of loops, and site 
perimeter 

In the same manner as in Section \Q we define the generating function counting 
the total value of the parameters j(A), £(A), p-(A) and p+(A) in animals on 
the triangular lattice A. Thus, let J7s(i), £s{t), Vg(t), Vtit) be the generating 
functions dealing with animals of source S, and ^T[s](t), C[s](t)i V^(t) an d 
■pjtjj(i) be their analogues counting animals with a source included in S. 

Animals on the triangular lattice A can be regarded as inflated heaps on 
the model (Q,C) defined in Section [U As on the square lattice, this leads us to 
define the following generating functions: 

• £s(t) counts heaps of base S marked with a piece at a position q such that 
either q — I or q + I is not in Q; 

• Us(t) counts heaps of base S marked with a piece at a position q such 
that q + 2 is not in Q; 

• if q — 2 is in S but q is not, W|(t) counts heaps of base S U {q}; 

• W S (t) is the sum of Wf (*) for all q in (5 + 2) \ S. 

Finally, define the counterparts of these generating functions that count heaps of 
base included in S, with Wr^, (t) counting heaps with a base included in 5U {q}, 
containing at least q — 2 and q. 

Again, the results of Section [3] enable us to compute all of these generat- 
ing functions. Thus, with the theorem below, one may compute the desired 
generating functions. 

Theorem 5.2.1. The generating functions counting the total number of adjacent 
sites, number of loops, and the lower and upper bound on the site perimeter in 
animals with a source included in S are: 

m' [s] (t)-iu [s] {t)-w [s] {t) 

J[s]{t) = — ~ t ; 

C [s] (t)=tJ [s] (t); (32) 

V^ s] (t) = \S\H [s] (t) + (t + t 2 )H' [s] (t) - J [s] (t); (33) 

(t) = 2\S\H [s] (t) + (t + t 2 )H{ s] (t) - £ [s] (t). (34) 

The generating functions counting animals with source S are: 

Js {t) = t 2 K's(t) - Msit) + 3{S)H s (t) - W s (t) _ 

1 ~\~ i 

C s (t) = tj s (t); (36) 
V s (t) - \S\H s (t) + (t + t 2 )H' s (t) - J s (t); (37) 
Vt<f) = 2\S\H s (t) + (t + t 2 )H' s (t) - C s (t). (38) 

First, let us show the identities {!!, (05), ([57]) . ([55 ]l . Let C [S ](t) (resp. C s (t)) 
be the generating functions counting the total number of sites only supported at 



29 



the center in animals with a source included in S (resp. of source S). Equation 
(ff2")) . found in Section implies that: 

C [s] (t)=t 2 H' ls] (t); 
C s (t)=t 2 H' s (t). 

Therefore, equations (IM)) , ([57)1 and (j3"5|) are found by summing the first 

identity of Lemma \b . 1 . 21 and idendity (|30p on all the relevant animals. 

Moreover, the external and internal site perimeters, like on the square lattice, 
are linked by: 

Vh(t) = V e s (t)-£ s (t). 

The case of animals with a source included in S is identical. 
The following lemma proves identities ([32| and |36|) . 

Lemma 5.2.2. The generating functions counting the total number of adjacent 
sites and loops are linked by: 

C [s] (t)=tJ [s] (t); 
C s {t)=tJ s {t). 

Proof. The proof is also similar to that of Lemma 14.2.21 for the square lattice. 
Let A be an animal of source S with two adjacent sites marked. We regard A 
as an inflated heap of the model (Q,C) marked with two adjacent pieces, and 
turn it into a factorized heap (Hi ■ H2) by pulling the marked pieces downwards, 
such that Hi has two maximal, adjacent pieces. 

We do the same with animals marked with a loop. The resulting factorized 
heaps differ only by a piece on top of Hi, which does not alter the base or 
neighbourhood of Hi. Therefore, Lemma [2.5.31 provides the desired link on the 
generating functions. □ 

To deal with the number of adjacent sites, we use the same intermediate as 
for the square lattice. Let Xp S j (t) (resp. I[ 3 s j) be the generating function counting 
factorized heaps (Hi ■ H2), such that H\ has two maximal pieces, with positions 
q and q + 2 (resp. q and q + 3) for some q. These functions are linked to J7[s] (t) 
in the following manner. 

Lemma 5.2.3. The following identity holds: 

J[s]{t) = TT^ W _ TTi^®' 

Proof. As the generating functions counting heaps and inflated heaps are the 
same, we may regard the factorized heaps counted by 2? s -,(t) and 2|g, (t) as 
factorized inflated heaps, i.e. factorized pre-heaps (Hi ■ H2) where Hi and H2 
are inflated heaps. 

Let (Hi- H2) be an inflated factorized heap with base included in S such that 
Hi has two maximal pieces x and y, which are at positions at distance 2 but 
are not adjacent (say, y is higher than x). We transform it by removing from 
Hi all pieces of the stack of y that are higher than x. The resulting factorized 
heap (H[ ■ H2) is such that H' x has two maximal pieces, either adjacent or at 
positions at distance 3 (see Figure ITB")). 
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Figure 13: Removing the pieces of the stack of y that are higher than x 
uncovers a maximal piece z, either at a position at distance 3 from x or 
adjacent to i. 



We now consider the inverse transformation. Let (Hi • H2) be an inflated 
heap with base included in S such that Hi has two maximal pieces at positions 
at distance 3. As S is aligned, Hi is also aligned, so its maximal pieces cannot 
be at the same height: let z be the higher piece and x the lower. The inverse 
transformation is performed by adding a stack of arbitrary height on top of z, 
at a distance 2 from x. 

Finally, let (Hi • H2) be a inflated heap of base included in S, such that 
Hi has two maximal, adjacent pieces. This time, the inverse transformation is 
performed by adding a stack of arbitrary height on top of either maximal piece, 
so that these heaps must be counted twice. We thus have: 



Zfs](*) 



J[ S ](t) = — t {Z[S](t) + 2J[S](t)), 



which yields the announced formula for i7rgi(t). 



□ 



We complete the proof of Theorem 15.2.11 by using the following result in 
conjunction with the results of Section [U 



Lemma 5.2.4. The following links hold: 



t 



1-t J 1 
t 



l-t V 



m[ s] (t)-u [s] (t) = M( s] (^-^ 

W[s\(t) = W [a] (j^- 



t 

Proof. The fourth identity is an instance of {2} (Section 12. 2| : it is due to the 
fact that every heap is obtained by replacing each piece of a strict heap by a 
stack of pieces. 

To prove the first three identities, we regard the heaps counted by U[s](t) 
as factorizcd heaps by pulling the marked piece downwards. In this way, these 
identities link generating functions counting factorized heaps and almost-stict 
factorized heaps. This link is found by replacing each piece of both parts of an 
almost-strict factorized heap by a stack of pieces. □ 
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Thus, combining identity ([25]) , Lemmas 14.3.21 and 14.3.31 and this result, we 
find: 



T 2 



(t) = -^— t (lf s] (t) + tH' [s] (t) - U [s] (t) - W [s] (*)) ; 
(-l) |SN - T| W [r] (t) = Ws(t) - j(S)Hs(t). 



TCS 

Lemma \5 . 2 . 31 and an inclusion-exclusion argument yield ([31]) and (|35|) . 

6 Applications 

6.1 Asymptotic results 

We now use Theorems 14.1.21 and 15.2.11 to compute asymptotic estimates of the 
average value of the various parameters in animals of a given area. Let S be a 
one-line source, and let /i be the growth constant of the animals of source S on 
the square lattice T. Propositions 13 . 2 .41 and [3~3~H1 assert that this constant does 
not depend on S, and that the growth constant of animals of source S on the 
associated triangular lattice A is p = p + 1. 

Moreover, let j(n), l(n), Pi(n) and p e (n) be the average number of adjacent 
sites, number of loops, internal and external site perimeter of animals of source 
S and area n on T; let j(n), l{n), p e (n), pi(n), p-{n) and p+(n) be the average 
value of each parameter of the animals on A. 

Corollary 6.1.1. Assume that T is cither the full lattice, the half lattice, or 
a cylindrical bounded lattice. As n tends to infinity, we have the following 
estimates: 

j{n) 7—; j(n) 



M+l M+l 

*(n) ~ 4; i(n) ~ — 

PiW ~ Pe(n) ~ p_ (n) ~ (l + ) n; 

1 



p+(n) ~ 1 



/2 + 1, 

Moreover, the average site perimeters on the triangular lattice verify for large 
n: 

p-(n) <pi(n) ~ p e (n) <p + {n). 

In the unbounded lattices, the growth constants are fi = 3 and p = 4. Thus, 
these estimates become: 

iOO ~ ^; pW ~ — ; JW ~ ^; ^-^^"20"' 

£(n)~- ; £( n )^_ ; p+ ( n )^_. 
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Proof. We use singularity analysis on the identities of Theorems 14. 1.2l and l5. 2. II 
Let us begin with the square lattice T. Let (Q,C) be the corresponding model 
of heaps of dominoes. As T is not a rectangular lattice, the position q + 2 is in 
Q as soon as q is. Thus, the generating function Us(t) is zero. 

According to Propositions 13.2.41 and I3.3.6[ the asymptotic behaviour of the 
coefficients of Hs(t) and Ws(t) differ only by a multiplicative constant; thus, 
the only meaningful term in Js(i) is Y^jH' s (t). This term yields the announced 
estimate of and the other parameters follow. 

All that is left to prove is that the term Es(t), needed to compute the total 
internal site perimeter, is negligible. If T is the full lattice or a cylindrical lattice, 
this term is zero. If T is the half-lattice, it counts heaps of base S marked with 
a piece at position 0: 

E s (t)=H { °\t). 
According to Proposition 1 2. 6. 31 we have: 

E s (t) < ^H [s] (t)H {0} (t), 

where the inequality holds coefficient-by-coefficicnt. Using Proposition l3.3.6l we 
see that this is equivalent to X^/l — 3t as t tends to 1/3, which yields negligible 
coefficients in comparison to H' s {t). 

The case of the triangular lattice is similar. □ 

In this corollary, the rectangular lattices are not included. Computing the 
estimates in this case is more difficult, because these lattices are the only ones 
where the generating functions Us it) and Es(t) have non- negligible coefficients 
compared to those of H' s {t). 

6.2 Examples 

6.2.1 Single-source animals on the full lattices 

In all the cases covered by the results of Sections 4 and 5, the simplest is the 
one of single-source directed animals on the square lattice. 

Corollary 6.2.1. The generating functions counting the total number of adjacent 
sites, number of loops and site perimeter of the single-source directed animals 
on the full square lattice are respectively given by: 

_1 / 1-4^ + 4*3 



2t(l + t)V VTTt(l - 3t) 3 / 2 
1> l-W + 4*3 



2 V V / TTI(1 - 3t)V 

1 ( „ . . . 2 . 1 - 3t + 2t 2 + 1 3 - 3t 4 



p (*) = w , ,A ~ 1 + t + t + 



2t(l+t)\ V / TTI(1 - 3t) 3 / 2 

The value of P(t) was conjectured by Conway [S], and the value of L(t) was 
proved by Bousquet-Melou using a gas model method [3]. 

Proof. To compute J(t), we use Theorem 14 . 1 . 2 1 with S — {0}. The generating 
function Usq\(£) is zero; W{o}(£) counts strict heaps of base {0,2}, that is, 
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animals with a compact source with two sites. The quantity j({0}) is zero. 
Lemma 13.3.31 thus entails that: 

_ t 2 A'(t)-D(t)A(t) 
J[) - 1 + t 

where the generating functions A(t) and D(t) are given by (|T7|) and (fTS")) . The 
values of L(t) and P(t) follow from Theorem [4.1.21 □ 

Corollary 6.2.2. The generating functions counting the total number of adjacent 
sites, numer of loops, and bounds on the site perimeter of single-source animals 
on the full triangular lattice are: 

1 / l-7t+ 12t 2 - 2t 3 
m = ^TT^(l-t (1 _4t)3/ 2 

1 -n+ m 2 - 2t 3 

(1 - At) 3 / 2 

2 1 - 6t + lit 2 - 52t 3 + 2t 4 

+ (1 - 4t)3/2 

3- llt + 8t 2s 



^=2^)(" 1 -* 



(1 - 4<)3/2 



This time, the value of C(t) is different from the one found by Bousquet- 
Melou, who used a different definition of loops. 



Proof. This time, we use Theorem 15.2. II Exactly as above, we have: 

t 2 A'(t) - V(t)A(t) 



Jit) 



l + t 



where A(t) and T>{t) are given by (fT8|) and (fl6|) . Once again, Theorem 15.2.11 
enables us to derive the three other generating functions. □ 

6.2.2 Compact-source animals on the full lattices 

According to Proposition 13. 3. 4"! the number of compact-source animals on both 
the square and triangular lattice is extremely simple. Below are the generating 
functions counting the total number of adjacent sites of these animals. 

Corollary 6.2.3. The generating functions counting the total number of adjacent 
sites of the compact-source directed animals on the full square and triangular 
lattices are given by: 

1 ( l-2t 1 - 3t - 2t 2 \ 

Jc(t) 



Jc(t) = 



2 \VTTt(l -3i)3/2 (l + t)(l-M) 2 

1 / l-3t 1 — 5* + 2t 2 \ 
2(1 +t) { (l-4<)3/2 ~ (l_4t) 2 1' 



Of course, although we only give the generating functions for the total num- 
ber of adjacent pieces, Theorems 14.1.21 and 1 5 . 2 . ll can still be used to deal with 
the total number of loops and total site perimeter. 
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Proof. Let k > 1; we use Theorem 14.1.21 on animals on the square lattice of 
source Ck (the compact source with k sites). The generating function Uc k {t) is 
zero. We remark that Wc k counts animals of source Ck+i; moreover, there are 
k — 1 pairs of adjacent sites in Ck- Thus: 

(t) - t2H ' c * {t) + (fc Z 1)gc * (t) Z g£*±iW 

On the triangular lattice, Theorem 1 5 . 2 . 1 1 entails that: 

t 2 H' Ck (t) + (k-l)Hc k (t)-H Ck+ M 



Jc h (t) = 



1 + t 



Lemma 13.3.31 gives the value of the generating functions Hc k (t) and Ttc k (t) ■ 
Summing these identities for all k > 1 and using the formula 

we obtain: 

Jc(t) = ^ (l?A' c {t) + {l A ^ [t)) 2 ~ A cH) - (A c (t) A(t))j . 

The generating function 3c® is computed in the same manner. Recall that 
the generating functions counting compact-source animals are: 

Ac® = ^ Ac® = U 

6.2.3 Single-source half-animals on the square rectangular lattices and half- 
lattice 

Let to > 1. We consider half-animals on the bounded rectangular lattice of 
width to (or animals with a single source at the corner of the lattice). Our 
goal is to compute the generating functions giving the average site perimeters 
of these half-animals, both internal and external. 

Corollary 6.2.4. The generating functions counting the total external and inter- 
nal site perimeter of half-animals in the rectangular square lattice of width to 
are: 

j£(t) - D m (t) + Y^D'm® + 7-^™w 2 ; 

PL® = j^- t D m (t) + ^—D' m ® - YT- t ( D m{t) - D m _ 2 (t))D m (t), 

where D m (t) is given in terms of the polynomials computed by Proposition l3 . 2 . D 

D ( t ) - ^Ziilfaj _ i 
lm \i+t) 



35 



Proof. The model of heaps of dominoes corresponding to the rectangular lattice 
of width m is Q = [0, m — 1]. Instead of studying pyramids of base in this 
model, we choose, by symmetry, to study pyramids of basis m — 1. These 
pyramids are still counted by D m (t), according to Proposition 13. 2. 21 

We now compute all the necessary generating functions to apply Theo- 
rem 13112] with S* = {to — 1}. As to + 1 is not in Q, the generating function 
W/ m _u is zero; the quantity j({m — 1}) is also zero. The two remaining gen- 
erating functions are given by: 

u {m . n (t) = 

where -ff[„_x}(*) counts strict heaps of base m — 1 marked with a piece at 
position g. 

As no heap marked with a piece can have an empty base, the generating 
functions H3_y.it) and H?r_y,At) are equal. We compute them using Propo- 
sition [HUSH for this, we need the generating functions counting strict heaps of 
base included in {to — 1} avoiding the positions m — 2 and 0. As no nonempty 
pyramid of base to — 1 can avoid to — 2, we have: 

Moreover, the pyramids of base m — 1 avoiding are the pyramids of the model 
Q = [2, to— 1], which is equivalent to the rectangular model with to— 2 positions. 
Thus, we get: 

V [{m-l}H 

The generating functions -ff{ TO _i}(£) and H{ y(t) are both D m (t). Finally, to 
compute Hi m -2}{t), we remark that a strict pyramid of base m — 1 is either 
reduced to a single piece at position m — 1, or is a piece at position m — 1 
surmounted by a strict pyramid of base to — 2. This entails that: 

H {m -i}{t) =t(l + H {m _ 2} (t)). 

Hence: 

tt m - Dm{t) i 

#{m-2}W — 7 I- 

Putting this all together, we get, using Proposition 12.6.31 
1 



Vr? m _ii 1 (t) = l + i?m-a(t). 



U{m-1}(*) 



1 + i 



T 



D m {t))D m (t) + D m (t) (^£^1 - i 



£{m-l}0) = — f(l + ^ m (t))i)mW+ (An(*)- D m ^{t))D m {t) 

_L — T" C- - 

The first identity simplifies into: 

D m {tf 



U{ m -l}(t) 

We conclude using Theorem 14. 1.21 □ 
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The value of the generating functions P^t) were conjectured by Le Borgne 
|10j . who gave an induction formula to compute them. By using the formula 
above, one can prove this conjecture. Finally, by simply letting m tend to 
infinity, this result extends to the unbounded half-lattice. 

Corollary 6.2.5. The generating functions counting the total external and inter- 
nal site perimeters of half-animals on the unbounded lattice are: 

P e (t) = D(t) + j^-D'it) + ^-D{t) 2 ; 
P\t) = ^ rt (D{t) + D'(t)). 
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